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This  investigation  considered  the  structure  and  mixing  properties  of  round  and 
plane  vertical  noncondensing  gas  jets  in  liquids,  motivated  by  application  involving 
metal  processing,  direct-contact  condensers,  gas  dissolution  systems,  reservoir 
distratification  systems,  liquid- metal  combustors,  and  nuclear  reactor  pressure- 
suppression  systems  -  among  others.  The  main  objective  was  to  provide  new 
measurements  of  flow  structure  and  mixing  properties,  however,  the  new 
measurements  were  also  used  to  evaluate  past  proposals  for  analyzing  the  process.  The 
first  phase  of  the  study  considered  initial  measurements  and  analysis  of  round  jets  (Loth 
&  Faeth,  1987);  the  present  investigation  involved  additional  measurements  and 
analysis  of  both  round  jets  and  plane  jets.  Supplementary  work  under  this  investigation 
considering  the  structure  of  gas  jets  in  gases,  has  also  been  reported  by  Chuech  (1987) 
and  Chuech  et  al.  (1987,  1988). 

Gas  jets  submerged  in  liquids  are  complicated  by  effects  of  unsteadiness  at  low 
flow  rates,  similar  to  the  pulsitile  behavior  of  highly  buoyant  single-phase  flows.  This 
involves  oscillatory  release  of  gas  which  can  cause  liquid  to  slug  into  the  jet  passage. 
The  unsteadiness  is  accompanied  by  appreciable  fluctuations  of  static  pressures  in  the 
flow  passage  and  beyond  the  jet  boundaries,  often  resulting  in  excessive  noise  and 
vibration  as  well  as  blockage  of  passage  in  cases  where  the  gas  reacts  with  the  liquid 
(Kemey  et  al.,  1972;  Chan,  1974;  Avery  &  Faeth,  1974). 

Effects  of  unsteadiness  can  be  reduced  by  increasing  gas  flow  rates  in  the 
passage,  frequently  leading  to  operation  with  underexpanded  jets,  where  the  flow  is 
sonic  and  the  static  pressure  is  greater  than  the  ambient  pressure  at  the  jet  exit.  It  is  well 
known  that  static  pressure  equalization  for  underexpanded  gas  jets  in  gases  occurs  in  an 
external-expansion  region  near  the  jet  exit  -  involving  shock  waves  and  other 
compressible-flow  phenomena  (Shapiro,  1954).  Moiseev  (1962)  and  Surin  et  al. 
(1983)  observed  that  dynamic  pressures  along  the  axis  of  underexpanded  gas  jets  in 
gases  and  liquids  were  similar  while  static  pressures  along  the  axis  of  underexpanded 
gas  jets  in  gases  and  liquids  were  also  found  to  be  similar  during  the  first  phase  of  this 
investigation  (Loth  and  Faeth,  1987).  These  observations  suggest  that  a  multiphase 
external-expansion  region,  involving  gas  dynamic  phenomena,  is  present  for 
underexpanded  gas  jets  in  liquids  as  well. 

In  view  of  the  complexities  of  unsteadiness  at  low  flow  rates,  and  the  presence 
of  the  external-expansion  region  at  high  flow  rates,  current  understanding  of  gas 
injection  into  liquids  is  not  very  complete;  nevertheless,  the  flow  has  received  some 


attention  in  the  past.  Earlier  work  in  this  laboratory  considered  round  jets,  including  the 
following:  noncondensing  gas  jets  in  liquids  (Tross,  1974);  condensing  gas  jets  in 
liquids  (Kemey  et  al.,  1972;  Weimer  et  al.,  1973;  Chen  &  Faeth,  1982);  and  reacting 
gas  jets  in  liquids  (Avery  &  Faeth,  1974;  Chen  &  Faeth,  1983).  The  earliest 
experimental  studies  were  confined  to  gross  parameters  like  the  length  of  the  vapor-or- 
gas-containing  region  of  condensing  or  reacting  jets  (Kemey  et  al.,  1972;  Weimer  et 
al.,  1973;  Avery  &  Faeth,  1974).  Analysis  was  also  undertaken  based  on  an  integral 
model  of  turbulence  and  the  locally-homogeneous  flow  (LHF)  approximation  of 
multiphase  flow  theory  (Soo,  1967).  The  LHF  approximation  implies  infinitely-fast 
interphase  transport  rates  and  local  thermodynamic  equilibrium  between  the  phases, 
i.e.,  velocity  differences  between  the  phases  are  neglected  and  phase  and  chemical 
equilibrium  are  assumed  to  be  maintained  at  each  point  within  the  flow.  Gas-dynamic 
processes  in  multiphase  flows  are  not  well  understood;  therefore,  the  complexities  of 
the  external-expansion  region  were  avoided  by  defining  effective-adapted-jet  exit 
conditions,  where  the  static  pressure  at  the  passage  exit  was  the  same  as  the  ambient 
pressure.  This  approach  was  reasonably  successful  for  correlating  the  length  of  the 
vapor-or  gas-containing  region  of  condensing  jets  (Kemey  et  al.,  1972;  Weimer  et  al., 
1973;  Avery  &  Faeth,  1974). 

Tross  (1974)  measured  mean  void  fractions,  using  an  electrical  conductivity 
probe,  and  dynamic  pressures,  using  an  impact  probe,  in  adapted  and  underexpanded 
air  jets  in  still  water.  The  flow  was  analyzed  similar  to  Weimer  et  al.  (1973)  and  Avery 
&  Faeth  (1974).  Unfortunately,  the  analysis  was  not  very  effective  for  predicting  flow 
structure  in  spite  of  its  earlier  success.  This  was  attributed  to  the  uncertainty  of  probe 
measurements  in  multiphase  flows  as  well  as  problems  of  treating  a  rapidly  developing 
flow,  having  a  large  density  ratio,  with  an  integral  model  (Tross,  1974). 

Subsequent  work  in  this  laboratory  considered  a  higher-order  turbulence  model 
in  an  effort  to  improve  predictions  (Chen  &  Faeth,  1982,  1983).  The  conserved-scalar 
formalism  of  Bilger  (1976)  and  Lockwood  &  Naguib  (1975)  was  used,  which  had 
provided  a  successful  treatment  of  constant  density,  variable  density,  and  combusting 
turbulent  jets  during  other  work  in  this  laboratory  (Jeng  &  Faeth,  1984).  The  LHF 
approximation  was  adopted,  while  only  flow  rates  slightly  in  excess  of  the  sonic 
condition  were  considered  to  avoid  problems  of  external  expansion.  A  key  feature  of 
the  combined  use  of  the  LHF  and  the  conserved-scalai  approximations  is  that  scalar 
properties  (void  fraction,  density,  etc.)  are  only  functions  of  mixture  fraction  (the 
fraction  of  mass  at  a  point  which  originated  from  the  injector).  These  functions,  called 
state  relationships,  can  be  found  from  straight  forward  adiabatic  mixing  or  adiabatic 
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reaction  calculations,  typical  of  fundamental  thermodynamic  theory.  Use  of  the 
turbulence  model  and  LHF  approximations  was  successful  for  predicting  the  length  of 
the  vapor-  or  gas-containing  region  of  both  condensing  and  reacting  jets,  with  all 
empirical  parameters  of  the  turbulence  model  fixed  at  values  appropriate  for  constant- 
density  single-phase  jets  (Chen  &  Faeth,  1982,  1983).  The  same  approach  has  been 
evaluated  for  a  variety  of  other  multiphase  jets,  see  Faeth  (1987)  and  references  sited 
therein.  In  general,  the  approach  was  reasonably  successful  in  dense  multiphase  flows 
but  tended  to  overestimate  the  rate  of  development  of  dilute  dispersed  flows  v/here 
finite  interphase  transport  rates  become  important  and  the  LHF  approximation  fails. 

Based  on  these  findings  the  same  approach  was  examined  during  the  first  phase 
of  the  present  investigation  (Loth  &  Faeth,  1987).  These  measurements  involved 
adapted  and  underexpanded  round  air  jets  injected  vertically  upward  in  still  water,  with 
underexpansion  ratios  (the  ratio  of  jet  exit  to  ambient  static  pressures  or  the  ratio  of 
initial  jet  mass  flow  rates  to  the  initial  jet  mass  flow  rate  for  sonic-adapted  conditions  at 
the  jet  exit)  up  to  8.  It  was  judged  to  be  premature  to  consider  the  complexities  of  a 
shock-wave-containing  multiphase  external  expansion  region;  therefore,  effective- 
adapted-jet  exit  conditions  were  used,  following  Kemey  et  al.  (1972).  This  was 
justified  by  recent  measurements  of  Birch  et  al.  (1984,  1987),  Chuech  (1987)  and 
Chuech  et  al.  (1987,  1988)  which  showed  that  effective-adapted-jet  exit  conditions 
provide  an  effective  treatment  of  the  structure  of  underexpanded  gas  jets  in  gases,  even 
quite  close  to  the  external  expansion  region  itself.  Because  of  the  rapid  development  of 
the  flow  for  gas  injection  into  liquids,  resulting  from  the  large  density  ratio  (Chen  & 
Faeth,  1982,  1983)  the  extent  of  the  external  expansion  region  should  be  smaller  for 
gas/liquid  than  gas/gas  systems,  favoring  the  approximation  as  well.  The  comparison 
between  predicted  and  measured  time-averaged  void  fraction  distributions  was 
moderately  successful  -  particularly  at  high  underexpansion  ratios.  For  adapted  and 
subsonic  exit  conditions,  however,  large-scale  and  pulsitile  jet  flow  was  observed 
which  enhanced  mixing  rates  well  above  predictions.  Furthermore,  unsteadiness  near 
the  jet  exit  appeared  to  influence  distributions  of  void  fractions  even  for  the  highest 
underexpansion  ratios  considered.  These  results  highlighted  the  need  to  study  flow 
properties  near  the  jet  exit  in  greater  detail. 

Related  studies  by  others  paralleled  activities  in  this  laboratory.  To  begin  with 
round  jets,  early  work  by  Cumo  et  al.  (1978),  Kudo  et  al.  (1974),  Lee  et  al.  (1979)  and 
Young  et  al.  (1974)  involved  measurements  of  only  the  gross  features  of  condensing 
jets,  like  the  length  of  the  vapor-containing  region.  Relevant  portions  of  this  data  were 
used  to  evaluate  LHF  analysis  of  condensing  jets  (Chen  &  Faeth,  1982).  Bakaklevskii 
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&  Chekhovich  (1978)  report  temperature  and  dynamic  pressure  profiles  in  condensing 
plane  jets  in  coflow,  however,  measurement  accuracy  and  flow  conditions  are  too 
uncertain  to  provide  more  than  a  rough  guide  to  structure.  Chan  (1974),  Chun  &  Sonin 
(1984),  Lambier  &  Chow  (1984)  and  Simpson  &  Chan  (1982)  measured  static 
pressures  near  the  passage  exit  for  subsonic  condensing  jets,  finding  the  large  pressure 
pulsations  indicative  of  unsteady  flow  that  were  mentioned  earlier.  Structure 
measurements  are  more  numerous  for  noncondensing  jets  but  most  results  are  limited  to 
dilute  bubbly  flows  far  from  the  passage  exit  (Abdel-Aal  et  al.,  1966;  Mahalingen  et  al., 
1976;  Ohba  et  al.,  1977  and  Ohba,  1979).  Past  studies  of  processes  near  the  passage 
exit  generally  involve  the  use  of  probes  or  other  intrusive  techniques,  for  example:  the 
dynamic  pressure  measurements  of  Moiseev  (1962)  and  Surin  et  al.  (1983); 
measurements  of  void  diameter  using  a  passage  plate  by  Bell  et  al.  (1972);  and 
measurements  of  entrainment  using  the  water  overflow  technique  by  Bell  et  al.  (1972) 
and  Carreau  et  al.  (1983). 

The  main  objective  of  the  present  investigation  was  to  extend  the  measurements 
of  Loth  &  Faeth  (1987)  for  round  jets  to  consider  the  near-jet-exit  region  of  plane  jets. 
Adapted  and  underexpanded  vertical  plane  jets  in  still  water  were  considered  with 
underexpansion  ratios  up  to  4.  The  plane  jet  configuration  offered  a  better  opportunity 
to  observe  the  flow  and  direct  photographs  and  shadowgraphs  were  obtained  for 
underexpanded  conditions.  Other  measurements  were  the  same  as  Loth  &  Faeth  (1987): 
time-averaged  void  fractions,  using  gamma-ray  absorption;  mean  and  fluctuating  jet  exit 
velocities  and  entrainment  rates,  using  laser-Doppler  anemometry;  and  static  pressures 
along  the  plane  of  symmetry  using  a  static-pressure  probe.  In  addition  hydrophone 
measurements  were  made  in  the  surrounding  liquid  near  the  jet  exit  to  study  acoustic 
phenomena  associated  with  unsteadiness  near  the  passage  exit.  The  present  structure 
measurements,  supplemented  by  earlier  measurements  of  Bell  et  al.  (1972),  Carreau  et 
al.  (1985)  and  Tross  (1974),  were  used  to  evaluate  predictions  based  on  past  proposals 
for  analyzing  this  flow,  e.g.,  use  of  the  LHF  and  effective-adapted-jet  approximations 
(Chen  &  Faeth,  1982,  1983;  Faeth,  1987;  Loth  &  Faeth,  1987). 

The  report  begins  with  a  description  of  experimental  and  theoretical  methods. 
Results  are  then  considered,  treating  flow  definition,  structure,  entrainment  and  the 
comparison  between  predictions  and  measurements  for  both  round  and  plane  jets.  The 
following  discussion  is  brief,  more  details  and  a  complete  tabulation  of  data  can  be 
found  in  Loth  (1988).  Additional  publications  relating  to  this  research  include:  Chuech 
(1987),  Chuech  et  al.  (1987,  1988)  and  Loth  &  Faeth  (1987,  1988). 


2.1 


Measurements  were  carried  out  in  a  water-filled  tank  (lm  x  2m  x  2m  high) 
which  had  glass  side  walls.  Figure  1  is  a  sketch  of  the  arrangement  for  tests  with  a 
round  jet.  The  jets  were  directed  vertically  upward  and  could  be  traversed  in  three 
directions  with  horizontal  and  vertical  positioning  accuracies  of  0.1mm  and  1mm.  The 
round  injectors  were  located  near  the  center  of  the  tank  and  were  generally  at  least 
650mm  below  the  liquid  surface  with  a  liquid  depth  of  1500mm.  In  this  case  flow 
disturbances  at  the  liquid  surface  were  controlled  by  wave  dampers  which  yielded  wave 
heights  less  than  75mm  and  pressure  disturbances  near  the  jet  exit  less  than  2kPa.  The 
plane  injector  was  located  near  the  surface  with  a  liquid  depth  of  1 100mm  and  wave 
dampers  were  not  needed. 

A  sketch  of  a  typical  round  injector  appears  in  Fig.  2.The  injector  consisted  of  a 
plenum,  a  honeycomb  flow  straightener  (1.6  mm  cells,  25  mm  long)  and  a  calming 
section,  all  having  a  diameter  of  32  mm,  followed  by  a  constant-radius  converging 
section  to  a  constant-diameter  flow  passage.  Two  jet  passages  were  used  having  inside 
diameters  of  4.9  and  1 1.0  mm  and  lengths  of  230  mm  to  yield  nearly  fully-developed 
pipe  flow  at  the  exit.  The  passages  extended  into  the  water  as  tubes  for  most  of  their 
length  and  were  tapered  from  the  outside  to  a  sharp  edge  at  the  exit.  Static  pressures 
were  measured  in  the  plenum  and  one  diameter  upstream  of  the  exit 

A  sketch  of  the  plane  injector  appears  in  Fig.  3.  This  consisted  of  a  140  cm3 
plenum  chamber  feeding  into  a  honeycomb  flow  straightener  (1.6  mm  cells,  20  mm 
long)  followed  by  three  30  mesh  brass  screens.  The  flow  then  passed  through  a  6:1 
lateral  contraction  designed  to  provide  a  uniform  exit  velocity  (within  1  percent  based 
on  potential  flow  theory).  The  dimensions  of  the  jet  exit  were  4.8  mm  wide  x  57  mm 
long,  yielding  an  initial  aspect  ratio  of  eleven.  Subsequently,  plate  glass  sidewalls 
provided  optical  access  for  a  streamwise  distance  of  13  exit  widths  and  a  lateral  distance 
of  16  exit  widths.  Screens  (30  mesh)  were  placed  along  the  sides  of  the  test  section  to 
reduce  disturbances  while  the  liquid  level  was  maintained  near  the  top  edge  of  the  test 
section. 

The  air  supply  was  drawn  from  laboratory  facilities  which  were  filtered  and 
dried  to  a  dewpoint  less  than  240  K.  Air  flow  rates  were  controlled  with  pressure 
regulators  and  metered  with  critical-flow  nozzles  having  diameters  of  2.43  and  7.62 
mm  (Flow  Engineering  Co.,  0.1  percent  accuracy).  Pressures  on  the  upstream  side  of 
the  critical-flow  nozzles  were  measured  with  a  250  mm  diameter  Heisse  pressure  gage 
(0-2  MPa  range). 


Figure  1.  Sketch  of  test  apparatus. 
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Figure  2.  Sketch  of  round  injector. 
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Figure  3.  Sketch  of  plane  injector. 
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2.2  Instrumentation 

Flow  Visualization.  Flash  photographs,  high-speed  motion  picture  photographs 
using  continuous  lighting,  and  flash  shadowgraphs  (plane  jet  only)  were  used  to 
document  the  mixing  and  dynamics  of  gas  release  at  the  jet  exit,  and  to  observe  the 
compressible  wave  field  for  the  plane  jets. 

Static  Pressures.  Static  pressures  were  measured  along  the  axis  of  the  1 1  mm 
diameter  round  jet,  and  along  the  center  plane  of  symmetry  of  the  plane  jets,  using  an 
arrangement  similar  to  Eggers  (1966).  The  arrangement  for  the  round  jet  is  illustrated  in 
Fig.  4. 

The  static  pressure  probe  consisted  of  a  1mm  diameter  stainless-steel  tube  with 
the  upstream  end  closed  and  a  0.4  mm  diameter  static-pressure  tap.  The  static-pressure 
tap  was  located  so  that  it  could  be  traversed  over  the  region  -  1  <  x/d  <,  7,  with  a 
positioning  accuracy  of  1  mm,  where  x  is  the  distance  from  the  jet  exit  and  d  is  the 
passage  exit  diameter.  The  probe  was  centered  along  the  axis  with  two  supports:  one 
within  the  passage  at  x/d  =  -  6,  the  other  in  the  bath  at  x/d  =  13.  A  small  air  purge  was 
used  to  keep  the  probe  free  of  water  when  measuring  air  jets  in  water.  The  arrangement 
for  the  plane  jet  was  similar  to  the  round  jet,  except  that  the  static  pressure  tube  was 
supported  upstream  of  the  converging  section  and  twenty-five  slot  widths  downstream 
of  the  exit  while  the  probe  could  be  traversed  in  the  range  -  1  <  x/b  <  13.  These  probes 
were  only  in  place  when  static  pressures  were  measured:  no  probes  were  located  in  the 
flow  field  for  all  other  operations.  Static  pressures  were  read  with  a  250  mm  diameter 
Heisse  pressure  gage  (0-2  MPa)  and  are  estimated  to  have  experimental  uncertainties 
(95  percent  confidence)  less  than  30  percent  of  the  pressure  difference  between  the  jet 
and  the  water  bath  (Loth,  1988). 

It  is  well  known  that  acoustic  feedback  can  cause  instabilities  which  enhance 
mixing  of  underexpanded  gas  jets  in  gases  (Sherman  et  al.,  1976).  This  affect  was 
studied  using  a  Spartan  Model  1 10-8459  Hydrophone  (2dB  response  in  the  range  1-40 
kHz)  during  tests  with  the  4.9  mm  diameter  round  passage  and  the  plane  passage.  The 
hydrophone  was  positioned  in  the  bath  at  x/d  =  6  and  r/d  or  y/b  =  6.  The  hydrophone 
output  was  processed  by  averaging  100  power  spectral  density  distributions  obtained 
with  a  rectangular  window  and  a  resolution  of  40  Hz. 

Void  Fractions.  Distributions  of  time-averaged  void  fractions  were  measured 
using  gamma-ray  absorption.  For  the  round  jets,  absorption  measurements  for  chord¬ 
like  paths  through  the  flow,  at  a  given  streamwise  distance,  were  deconvoluted  to  find 
radial  distributions  of  time-averaged  void  fraction.  Gamma-rays  were  obtained  from  a 
Co-57  source  (Amhersham  Inc.,  2  mCi,  271  day  half-life)  placed  in  a  lead  casket,  with 
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an  output  beam  of  gamma-rays  directed  across  the  flow  using  a  lead  collimator  (2  mm 
diameter  and  13  mm  long).  After  traversing  the  flow,  the  beam  passed  through  a  lead 
aperature  (2,  3,  and  5  mm  diameter,  depending  on  the  jet  width,  and  13  mm  long)  to 
the  detector  (EG  &  G  Model  905-1  sodium  iodide  scintillator  and  photomultiplier).  The 
detector  signal  was  preamplified  and  then  processed  by  an  EG  &  G  Model  590A  single¬ 
channel  analyzer  and  amplifier,  and  an  EG  &  G  Model  974  timer/counter.  The  output 
counts  of  the  timer/counter  were  then  collected  and  stored,  for  various  sampling  times, 
using  an  IBM-AT  computer.  The  absorption  window  of  the  detector  was  1 14-128  keV, 
to  capture  the  122  keV  emission  of  the  source  while  eliminating  low  energy  background 
radiation.  The  source  and  detector  were  located  on  either  side  of  the  jet  flow, 
submerged  in  the  tank  in  waterproof  housings.  At  each  axial  position,  absorption 
measurements  were  obtained  for  20-40  horizontal  chord-like  paths  through  the  flow. 
Roughly  10,000  counts  were  accumulated  for  each  path  which  required  sampling  times 
of  8-50  minutes.  The  linear  absorption  coefficient  of  gamma-rays  for  air  is  about  four 
orders  of  magnitude  smaller  than  for  water,  therefore,  the  logarithm  of  the  counting  rate 
is  proportional  to  the  fraction  of  air  in  the  radiation  path.  Assuming  axisymmetric  flow, 
the  absorption  measurements  were  deconvoluted,  following  Santoro  et  al.  (1981),  to 
obtain  radial  profiles  of  time-averaged  void  fraction. 

The  arrangement  for  the  plane  jet  measurements  was  very  similar  to  the  round 
jets,  except  that  the  source  was  replaced  by  a  10  mm  Ci  source  to  shorten  test  times. 
Naturally,  in  this  case  the  measurements  could  be  used  directly  and  no  deconvolution 
was  necessary. 

Bias  errors  of  the  void  fraction  measurements  are  influenced  by  whether  liquid 
laminae  in  the  absorption  path  are  parallel  or  normal  to  the  path  (Ohba,  1979).  Low 
intensity  ratios  reduced  uncertainties  due  to  this  effect  to  less  than  5  percent,  based  on 
the  findings  of  Ohba  (1979).  Dynamic  bias  errors,  due  to  turbulent  fluctuations  of  void 
fractions  along  the  path,  were  estimated  by  deconvolutions  through  a  linearized  error 
analysis,  which  showed  bias  errors  of  less  than  5  percent 

The  main  source  of  experimental  uncertainty  for  the  round  jets  was 
amplification  of  count  rate  uncertainty  during  deconvolution.  This  was  calculated  using 
the  matrix  model  of  Limbaugh  &  Kneile  (1984)  for  results  obtained  by  deconvolution 
of  data  with  a  prescribed  uncertainty.  This  resulted  in  experimental  uncertainties  (95 
percent  confidence)  at  the  centerline  of  less  than  20  percent,  with  relative  uncertainties 
increasing  inversely  proportional  to  the  indicated  mean  void  fraction  at  off-axis 
positions.  Since  no  deconvolution  was  necessary  for  the  plane  jet;,  their  experimental 
uncertainties  were  lower,  less  than  10  percent.  Assessment  of  errors,  based  on 


retrieving  a  mean  void  fraction  of  unity  for  several  positions  in  the  jet  air  core, 
confirmed  uncertainty  estimates  within  these  limits  (Loth,  1988). 


Velocities.  A  conventional  dual  beam,  frequency- shifted,  forward-scattev,  laser- 
Doppler  anemometer  (LDA)  was  used  to  measure  mean  and  fluctuating  streamwise 
velocities  at  the  injector  exit  (with  the  tank  empty)  and  mean  jet  entrainment  velocities  in 
the  liquid.  The  uncertainties  of  the  entrainment  velocity  measurements  (95  percent 
confidence)  are  estimated  to  be  less  than  10  percent  (Loth,  1988). 

Entrainment.  Jet  entrainment  rates,  defined  as  the  rate  increase  of  the  mass  flow 
rate  of  the  jet  with  the  streamwise  distance,  were  obtained  from  the  entrainment  velocity 
measurements.  These  measurements  were  made  along  the  jet  boundaries  outside  the 
gas-containing  region  at  a  radial  distance  r*,  or  lateral  distance  y*,  while  correcting  for 
low  levels  of  streamwise  velocity,  Um,  in  the  region  where  accurate  measurements  of 
radial  velocity,  Voo,  could  be  made.  Based  on  these  measurements,  the  rate  of  increase 
of  the  mass  flow  rate  of  the  jet,  m  with  streamwise  distance  is: 

dm/dx  =  (27tToo)jpoo(u«,  tan  J3oo  -  Voo);  j  =  1  (round)  =  0  (plane)  (2. 1) 

where  p  denotes  density,  poo  is  the  angle  between  the  jet  boundary  and  the  axis,  and  an 
overbar  denotes  a  time  average.  From  this  a  dimensionless  entrainment  coefficient  can 
be  defined  using  scaling  laws,  as  follows: 

Ce  =  dm/dx/  (Fe  poo  bH) 1/2  (2.2) 

where  Fe  is  the  streamwise  thrust  on  the  passage  of  the  round  jet,  and  is  the  thrust  per 
unit  length  for  the  plane  jet 

The  term  involving  Uoo  in  equation  (2.1)  contributed  less  than  10  percent  to  the 
entrainment  while  entrainment  rates  were  relatively  independent  of  the  location  of  y,* 
outside  the  gas-containing  region.  Estimated  experimental  uncertainties  (95  percent 
confidence)  of  the  entrainment  measures  were  less  than  20  percent  for  the  round  jets 
and  less  than  10  percent  for  the  plane  jets,  largely  governed  by  uncertainties  of  the 
velocity  measurements  and  effects  of  locating  r«,  or  y«>  (Loth,  1988). 

2.3  Test  Conditions 

Test  conditions  are  summarized  in  Table  1  as  a  function  of  underexpansion 
ratio,  m/ms.  For  the  round  jets,  flow  rates  varied  from  subsonic  jets  to  highly 
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Table  1  Summary  of  Test  Conditions-4- 


pe(kg/m3) 

me(g/s) 

Re/105 1 

Ri/10-4  § 

4.9  mm  iniector:  ffullv-developed  pipe  flowl: 

0.6 

1.36 

5.28 

1.11 

0.80 

8.14 

1.0 

1.52 

8.68 

2.73 

1.34 

3.31 

2.0 

3.06 

17.5 

7.42 

3.62 

1.21 

4.0 

6.12 

35.8 

16.7 

8.01 

0.54 

8.0 

12.3 

71.7 

35.9 

17.2 

0.25 

ILQ.mm. 

iniector:  ifullv-developed  pipe  flow'): 

0.6 

1.36 

26.3 

5.53 

1.80 

18.5 

1.0 

1.52 

43.8 

13.8 

3.01 

7.41 

2.0 

3.06 

87.7 

37.3 

8.16 

2.73 

4.0 

6.12 

174.0 

83.6 

18.5 

1.22 

4.8  mm  plane  iniector;  (slug  flow): 

1.0 

1.45 

113 

35.5 

1.22 

3.25 

2.0 

2.90 

226 

96.3 

3.30 

0.88 

3.0 

4.36 

339 

157 

5.40 

0.50 

4.0 

5.82 

453 

218 

7.49 

0.34 

£ 


p* 


'♦'Initial  conditions  for  an  air  jet  injected  vertically  upward  in  still  water  having  a 
temperature  of  294  ±  2K;  ambient  pressure  at  jet  exit  of  107.5  ±  0.7  kPa.  Air 
stagnation  temperature  of  295  =  2  K. 

^e  =  peF ed/(^e  me). 

§Ri  =  (poo  /  pe  -  1)  ad(me  /  Fe)2. 
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underexpanded  sonic  jets  (0.6  <  m/ms  <  8.0).  Measurements  for  the  plane  jets  all 
involved  sonic  exit  conditions  (1  <  m/ms  <  4). 

For  m/ms  >  1,  the  ratios  of  the  passage  exit  to  ambient  static  pressures  are  equal 
to  m/ms.  Flow  Reynolds  numbers  are  quite  high,  ca.  105  and  the  jets  were  turbulent. 
Richardson  numbers  were  relatively  low,  ca.  10*4;  therefore,  effects  of  buoyancy  were 
small  in  the  external-expansion  region  near  the  passage  exit.  However,  mixing  is  rapid 
for  air  injection  into  water  and  the  large  density  difference  between  water  and  air  causes 
significant  effects  of  buoyancy  farther  downstream,  x/d  >  20,  for  the  round  jets. 

Related  measurements  from  earlier  work  were  also  used  to  examine  flow 
structure.  These  measurements  all  involved  round  jets  as  summarized  in  Table  2.  The 
test  conditions  fall  within  the  range  of  present  measurements  for  round  jets,  although 
measurements  were  generally  made  using  intrusive  methods. 
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'*.■  Table  2  Test  Conditions  for  Related  Measurements3 
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$ 

Source 

Bell  et  al. 

(1972) 

Carreau  et  al. 
(1985) 

Tross 

(1974) 

R 

Injected  Fluid 

Nitrogen 

Nitrogen 

Air 

■ 

Jet  Direction 

Horizontal 

Horizontal 

Vertical 

Jet  Diameter  (mm) 

0.5  -  1.7 

0.5 

3.2  &  6.4 

m/ms 

3.5 

2-50 

1.4  -  14.0 

Rex  105 

0.5  -  2.0 

0.3  -  70 

0.8  -  5.0 

V 

V 

Measurements 

Entrainment: 

Water  Overflow 

Entrainment: 

Water  Overflow 

Void  Fraction: 

Resistivity  Probe 

Jet  Width: 

Dynamic  Pressure: 

Ji 

Passage  Plate 

Impact  Probe. 

Analysis  was  similar  to  the  approach  proposed  by  Chen  &  Faeth  (1982,  1983) 
to  treat  flows  of  this  type.  The  method  has  been  evaluated  in  this  laboratory  using 
structure  measurements  for  a  wide  range  of  multiphase  round  jets  (see  Faeth  (1987)  and 
references  sited  therein):  the  objective  of  present  work  was  to  extend  this  evaluation  to 
underexpanded  gas  jets  in  liquids. 

The  three  major  assumptions  of  the  analysis  are  as  follows:  (1)  use  of  the 
locally-homogeneous  flow  approximation  to  treat  multiphase  flow  effects;  (2)  use  of 
effective-adapted-jet  exit  conditions  to  treat  the  external-expansion  region  of 
underexpanded  jets;  and  (3)  use  of  a  higher-order  turbulence  model  to  treat  turbulent 
mixing.  The  LHF  approximation  implies  negligible  relative  velocities  between  the 
phases  and  local  thermodynamic  equilibrium;  therefore,  the  flow  is  treated  like  a  single¬ 
phase  fluid  having  large  density  variations  due  to  changes  in  gas  concentrations  while 
separated-flow  parameters,  like  drop  and  bubble  size  distributions,  do  not  enter  the 
formulation.  Consistent  with  present  limited  knowledge  concerning  the  structure  of  gas 
jets  in  liquids,  it  seems  prudent  to  evaluate  the  performance  of  LHF  analysis,  as  a 
baseline,  before  undertaking  the  additional  complications  of  separated-flow  analysis. 
Furthermore,  recent  evaluations  of  the  LHF  approximation  suggests  reasonably  good 
performance  in  the  near-injector  region  of  round  sprays  and  bubbiy  jets  (Faeth,  1987; 
Ruff  et  al.,  1988;  and  Sun  &  Faeth,  1986):  this  provides  additional  motivation  for 
examination  of  the  approach  here. 

The  effective-adapted-jet  approximation  is  frequently  used  to  avoid  the 
complexities  of  treating  gas-dynamic  phenomena  in  external-expansion  regions  when 
estimating  turbulent  mixing  for  both  single-  and  multi-phase  flows  (Avery  &  Faeth, 
1974;  Birch  et  al.,  1984,  1987;  Chuech  etal.,  1988;  Kemeyetal.,  1972;  Weimeretal., 
1973).  Recent  evaluations  have  also  shown  that  the  approach  is  reasonably  effective  for 
estimating  the  structure  of  underexpanded  gas  jets  in  gases,  even  close  to  the  external- 
expansion  region  (Birch  et  al.,  1984,  1987;  Chuech,  1987,  Cheuch  et  al.,  1987,  1988). 
Since  present  measurements,  and  those  of  Surin  et  al.  (1983),  suggest  that  external- 
expansion  regions  of  underexpanded  gas  jets  in  liquids  are  smaller  than  for 
underexpanded  gas  jets  in  gases,  the  region  where  the  approximation  is  useful  should 
be  even  larger  for  the  present  flows. 

Due  to  the  high  Reynolds  numbers  of  present  flows,  some  degree  of  modeling 
must  be  accepted  to  treat  their  mixing  proper  ies.  Consistent  with  past  work  on  gas  jets 
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in  liquids  and  related  multiphase  jets  in  this  laboratory  (Chen  &  Faeth,  1982,  1983; 
Faeth,  1987;  Ruff  et  al.,  1988;  Sun  &  Faeth,  1986),  turbulent  mixing  was  treated  using 
a  simplified  k-£-g  turbulence  model.  The  approach  is  similar  to  an  early 
recommendation  by  Lockwood  &  Naguib  (1975)  but  is  extended  to  use  mass-weighted 
(Favre)  averages  to  simplify  treatment  of  density  fluctuations  as  recommended  by 
Bilger  (1976).  Empirical  constants  for  round  jets  were  based  on  calibrations  for 
constant-  and  variable-density  single-phase  round  jets  using  the  Favre-averaged 
formulation  (Jeng  &  Faeth,  1984),  however,  the  values  are  not  very  different  from 
Lockwood  &  Naguib  (1975).  Empirical  parameters  were  the  same  for  plane  jets,  except 
for  use  of  a  modified  turbulent  Prandtl/Schmidt  number  as  recommended  by  Lai  et  al. 
(1986).  This  approach  is  plausible  since  the  present  flows  are  geometrically  simple 
high  Reynolds  number  turbulent  flows  for  which  turbulence  models  were  developed 
and  have  been  reasonably  successful.  Furthermore,  recent  evaluation  of  this  approach 
for  the  near-injector  region  of  large-scale  pressure-atomized  sprays  has  been  reasonably 
successful  (Ruff  et  al.,  1988)  suggesting  some  capability  to  treat  large  density  ratios, 
comparable  to  the  present  flows,  as  well. 

Other  major  assumptions  of  the  analysis  are  as  follows:  steady  (in  the  mean) 
axisymmetric  turbulent  jet  w;th  no  swirl;  boundary-layer  approximations  apply;  equal 
exchange  coefficients  of  all  species,  phases  and  heat;  buoyancy  considered  in  the 
governing  equations  for  mean  quantities  but  ignored  in  the  governing  equations  for 
turbulence  qualities;  and  effects  of  mean  kinetic  energy,  viscous  dissipation  and 
compressibility  are  ignored.  The  first  assumption  is  a  condition  of  the  experiments 
while  use  of  the  boundary-layer  approximations  is  not  an  issue  for  nonswirling  jets. 
Assuming  equal  exchange  coefficients  of  all  species,  dispersed-phases  and  heat  is 
widely  accepted  for  high  Reynolds  number  turbulent  flows  where  turbulent  transport 
dominates  flow  properties  (Lockwood  &  Naguib,  1975;  Bilger,  1976;  Jeng  &  Faeth, 
1984;  Faeth,  1987).  Neglecting  buoyancy  in  the  governing  equations  for  turbulence 
quantities  minimizes  empiricism  while  past  work  shows  that  the  effect  of  this 
approximation  is  small  with  respect  to  mean  properties  (Jeng  &  Faeth,  1984).  Finally, 
underexpanded  jets  involve  supersonic  How  downstream  of  the  passage  exit  so  that 
neglecting  effects  of  mean  kinetic  energy,  viscous  dissipation  and  compressibility  is 
clearly  questionable.  However,  the  external  expansion  region  is  not  large  and  is  only 
treated  approximately  using  the  effective-adapted-jet  approach;  therefore,  the  added 
complications  of  treating  these  effects  are  not  justified. 
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3.2  Formulation 

Under  the  assumptions  of  the  analysis,  the  conserved-scalar  formalism  can  be 
used  (Bilger,  1976).  This  involves  solving  governing  equations  for  conservation  of 
mass,  mean  momentum,  mean  mixture  fraction,  turbulence  kinetic  energy,  the  rate  of 
dissipation  of  turbulence  kinetic  energy,  and  mixture  fraction  fluctuations  squared.  All 
quantities  are  formulated  in  terms  of  Favre  averages,  which  are  defined  as  follows: 

<i>  =  p4>/p  (3.1) 

where  an  overbar  indicates  a  conventional  time  average.  By  this  definition  <|>  =  <j>  =  <j>", 
where  p4>"  =  0,  but  <j>"  *  0  necessarily,  unless  the  density  is  constant. 

The  Favre-averaged  governing  equations  can  be  written  in  the  following  form: 

r)0/9x  (pu  <)>)  +  d/dr  (dpv0)  =  d/dr  ((dpt/Cty)  d<p/dr)  +  dS^  (3.2) 

where  j  =  1  for  round  jets  and  j  =  0  for  plane  jets.  While  lateral  variations  are  denoted 
3(  )/dr  for  the  plane  jet  by  this  notation,  they  should  be  recognized  to  imply  d(  )/dy.  The 
dependent  variable  <)>  is  a  generic  variable  where  <j>  =  1  (for  conservation  of  mass),  u,  f , 
k,  e  and  g  (where  k,  e  and  g  are  Favre-averaged  quantities  in  this  formulation).  The 
source  terms,  S<|>,  and  the  appropriate  empirical  constants  are  summarized  in  Table  3. 
As  noted  earlier,  the  only  difference  in  the  round  and  plane  jet  formulation  involves 
changing  Of  =  ag  from  0.7  (for  round  jets)  to  0.5  (for  plane  jets)  as  recommended  by 
Lai  et  al.  (1986).  The  turbulent  viscosity  is  calculated  in  the  usml  manner: 

Ht  *  C^pk2/e  (3.3) 
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Table  3.  Summary  of  Source  Terms  in  the  Governing  Equations 
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M-  +  (M-t/Cfe) 

[CeiPt  (d  u/dr)2  -  CE2  p  e]  e/k 
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(p/Sc)  +  (|it/ag) 
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Cji  Cei 

Cgi  Ce2  =  Cg2 
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0.09  1.44 

2.8  1.87 
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3.3  Initif 


Several  effective-adapted-jet  approximations  have  been  proposed  to  treat  the 
mixing  properties  of  underexpanded  jets  while  avoiding  the  complications  of  the 
external-expansion  region,  see  Chuech  et  al.  (1988)  for  a  review  of  methods  proposed 
thus  far.  The  divergent-nozzle  approach  of  Kemey  et  al.  (1972)  was  used  during  the 
present  investigation.  This  involves  replacing  the  actual  external-expansion  process  by 
isentropic  flow  to  the  ambient  pressure,  and  applying  the  new  diameter,  velocities,  etc., 
of  the  flow  at  the  exit  plane  of  the  actual  passage  -  ignoring  the  presence  of  any  virtual 
origin.  This  approach  has  been  effective  for  estimating  the  mixing  properties  for 
injection  of  gases  into  gases,  and  conserves  mass,  momentum  and  energy  for  the  flow 
(Chuech  et  al.,  1988).  The  character  of  the  flow,  either  fully-developed  pipe  flow  or  a 
low- turbulence  (slug)  flow,  was  preserved  in  this  process  as  recommended  by  Chuech 
et  al.  (1988),  however,  use  of  either  slug  flow  or  fully-developed  flow  had  a  negligible 
effect  on  computed  results  in  the  region  where  measurements  were  made.  Formulas  for 
these  computations  are  reported  by  Loth  &  Faeth  (1987)  and  are  also  widely  available 
in  the  literature,  e.g.,  see  Shapiro  (1954).  This  approach  will  be  denoted  the  equivalent 
exit  condition  (EEC)  model  in  the  following. 

By  definition  f  =  1  and  g  =  0  across  the  jet  exit,  yielding  standard  initial 
conditions.  Present  measurements,  however,  suggested  enhanced  mixing  very  close  to 
the  jet  exit  due  to  effects  of  intrinsic  unsteadiness.  To  represent  this  effect,  a  partially- 
mixed  initial  condition  was  also  used  to  obtain  a  better  fit  of  void  fraction  distributions 
and  entrainment  rates  near  the  jet  exit.  The  mixed  initial  condition  involved  a  mixing 
layer  having  a  sinusoidal  variation  of  time-averaged  void  fraction  at  the  outer  edge  of 
the  injected  fluid  as  sketched  in  Fig.  5.  The  dimensionless  thicknesses,  Sic.  were  38 
percent  and  28  percent  for  the  round  and  plane  jets.  The  time-averaged  void  fraction 
distributions  were  found  from  the  measurements,  assuming  similarity  of  the  growth 
rate  of  mixing  layer  thickness  and  extrapolating  to  x/d  or  b  of  0.25.  The  associated 
velocities,  mixture  fractions  and  densities  were  computed  based  on  conservation  of 
momentum  for  a  constant-pressure  mixing  layer  and  assuming  Dirac  delta  functions  at  f 
=  0  and  1  for  the  probability  density  function  (PDF)  of  mixture  fraction. 

Measurements  of  jet  exit  conditions  for  the  round  jets  showed  that  fully- 
developed  turbulent  pipe  flow  was  nearly  achieved  (Loth  &  Faeth,  1987).  Thus,  mean 
velocity  distributions  at  the  jet  exit  used  the  power  law  distributions  due  to  Nikuradese 
(Hinze,  1975),  while  measurements  due  to  Laufer  were  used  to  specify  k  and  £ 


(Schlichting,  1979).  The  plane  jet  initial  condition  approximated  slug  flow  and  was 
taken  to  be  a  uniform  velocity  with  a  turbulence  intensity  of  less  than  3  percent  for  99 
percent  of  the  flow  width. 


3.4  Stat  e  Relationships 

When  exchange  coefficients  are  equal,  and  effects  of  kinetic  energy  and 
radiation  are  small,  all  instantaneous  scalar  properties  (density,  temperature,  phase 
fractions,  species  concentrations,  etc.)  are  only  a  function  of  mixture  fraction.  This 
implies  that  instantaneous  scalar  properties  can  be  found  by  straightforward  adiabatic 
mixing  (or  chemical  equilibrium)  calculations,  where  f  kg  of  passage  exit  fluid  and  (1-f) 
kg  of  ambient  fluid  are  adiabatically  mixed  and  brought  to  thermodynamic  equilibrium. 
The  relationships  between  scalar  properties  and  f  are  termed  state  relationships:  several 
examples  of  state  relationships  and  their  construction  appear  in  the  literature  (Bilger, 
1976;  Chen  &  Faeth,  1982, 1983;  Faeth,  1987;  Jeng  &  Faeth,  1984). 

In  the  present  case,  state  relationships  constructed  for  isothermal  mixing  of  air 
and  water  neglecting  the  small  vapor  pressure  of  water  were  found  to  be  adequate; 
therefore,  only  this  formulation  will  be  shown.  The  mass  fractions  of  air  and  water  are 
given  as 


Ya  =  f,  Yw  =  1-f  (3.4) 

For  isothermal  mixing,  the  phase  densities  remain  constant,  therefore,  the  mixture 
density  becomes: 


p  =  (f/pa  +  (l-f)/pw)’1 


(3.5) 


Finally,  the  void  fraction  is 


a  =  pf/pa  (3.6) 

In  order  to  evaluate  potential  effects  of  additional  property  variations  at  high 
underexpansion  ratios,  state  relationships  were  also  constructed  assuming  mixing  of  air 
at  the  maximum  effective  jet  exit  condition  considered  during  this  investigation  m/ms  = 
8).  Effects  of  kinetic  energy  were  also  ignored  for  this  calculation,  with  the  static  state 
at  the  passage  exit  taken  to  be  equal  to  the  ambient  pressure.  This  results  in  low  mixture 
temperatures  so  that  the  formation  of  ice  must  be  considered. 


State  relationships  for  the  isothermal  mixing  and  maximum  isentropic  expansion 
states  are  illustrated  in  Fig.  6.  Differences  between  the  two  state  relationships  are  only 
evident  at  mixture  fractions  near  unity,  where  ice  is  formed  for  isentropic  expansion,  as 
noted  earlier.  During  mixing  calculations,  flow  conditions  dropped  rapidly  below 
mixture  fractions  where  the  two  state  relationships  differ;  therefore,  the  simple 
isothermal  mixing  approximation  was  adopted  with  little  error  (Loth  &  Faeth,  1987). 

The  results  illustrated  in  Fig.  6  show  that  a  is  a  very  nonlinear  function  of 
mixture  fraction,  increasing  rapidly  with  increasing  f  near  f  =  0  and  then  remaining  near 
unity  thereafter.  This  behavior  can  be  quantified  by  noting  that  equation  (3.6) 
approaches  the  following  form  for  small  f: 

a  =  fpw/pa  ,  for  f  «  pa/pw  (3.7) 

Since  r>w/pa  ~  103  for  present  conditions,  equation  (3.7)  implies  that  a  is  a  very 
sensitive  function  of  f  in  this  region. 

3.5  Scalar  Properties 

Given  state  relationships,  time-averaged  and  Favre-averaged  scalar  properties 
were  computed  according  to  Bilger  (1976).  This  involves  finding  a  Favre-averaged 
probability  density  function  of  mixture  fraction,  P(f),  as  described  later.  Given  P(f),  the 
Favre-averaged  mean  and  mean-squared  fluctuating  values  of  any  scalar  property 
become  (Bilger,  1976): 


p  =  J<l>(f)P(f)df 

(3.7) 

$"2  =  /(<)>  (f)-$)2  P(f)  df 

(3.8) 

where  <j)(f)  is  the  state  relationship  for  the  property  <J>.  The  time-averaged  probability 
density  function  of  mixture  fraction,  P(f),  is  related  to  P(f),  as  follows: 


P(f)  =  p  P(f)  /  P(f) 


(3.9) 


Thus,  time-averaged  mean  and  mean-squared  fluctuating  values  of  any  scalar  property 
become: 


<t>  =  /  <l>  (0  P(f)  df  =  p  I  (p(f)  /  p(f))  P(f)df 


(3.10) 
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4>'2  =  P  I  ((<K0 - <t>)2 /  P(f))  P(f)  df  (3.11) 

The  time-averaged  density,  p,  which  is  needed  to  solve  the  governing  equations,  is 
found  by  setting  <J>  =  1  in  equation  (3.10). 

To  complete  the  formulation,  a  functional  form  must  be  assumed  for  P(f), 
although  the  specific  form  does  not  have  a  strong  effect  on  predictions  (Lockwood  & 
Naguib,  1975).  Any  two  parameter  distribution  can  be  used,  within  the  present  level  of 
closure,  since  its  two  unknown  parameters  can  be  found  by  noting: 

f  =1  f  P(f)  df 
g  =  !  (f-f)2  P(f)df 

Since  f  and  g  are  known  from  the  solution  of  the  governing  equations,  equations 
(3.12)  and  (3.13)  provide  two  implicit  expressions  to  solve  for  the  two  parameters  of 
the  PDF. 

Formally,  the  state  relationships  of  equations  (3.4)-(3.7)  imply  that  air  and 
water  never  coexist  at  a  point:  either  pure  water  or  pure  air  is  present.  Thus  the  correct 
state  relationship  for  the  flow  involves  Dirac  delta  functions  at  f=0  and  1,  which  were 
used  in  specifying  initial  conditions  as  noted  earlier.  As  a  practical  matter,  however, 
spatial  resolution  of  all  measurements  considered  during  this  study  was  not  sufficient  to 
resolve  small  dispersed-phase  elements  like  drops  and  small  bubbles;  therefore,  it  is 
also  realistic  to  consider  heterogeneous  mixtures  of  both  phases  in  any  realizable 
measuring  volume  for  observations.  Thus,  a  clipped-Gaussian  probability  density 
function  was  used  during  the  present  investigation,  following  Lockwood  &  Naguib 
(1975)  and  earlier  work  in  this  laboratory  (Chen  &  Faeth,  1982,  1983;  Jeng  &  Faeth, 
1984).  Solution  of  the  resulting  transcendental  equations  for  the  most  probable  value 
and  variance  of  this  distribution  was  facilitated  by  a  table  constructed  by  Shearer  et  al. 


(3.12) 

(3.13) 
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(1979). 

The  formulation  for  scalar  properties  is  actually  simpler  for  the  more  correct 
version  where  the  PDF  of  mixture  fraction  consists  of  Dirac  delta  functions  at  f  =  0  and 
1.  This  approach  will  be  stated  for  completeness  although  it  is  only  suitable  for 
noncondensing  and  nonreacting  flows.  Denoting  scalar  properties  at  the  jet  exit,  where 
f  =  1,  to  be  <j>0;  and  scalar  properties  in  the  ambient  environment  of  the  jet,  where  f  =  0, 


to  be  4>c«;  the  various  Favre-and  time-averaged  scalar  properties  using  the  Dirac  delta 
function  PDF  become: 


P  =  poo  (1  -  f  )  +  4>0  f 

(3.14) 

P  =  (poo  po  (1  -  f)  +  Po  poo  f)  /  (Pod  -  f)  +  poo  f) 

(3.15) 

P"2  =  (Poo-p0)2f  (I'd 

(3.16) 

p12  =  (poo  -  Po)2  P2  f  (1  -  f )  /  (poo  Po) 

(3.17) 

where 


P  =P0  Poo  /  (p0  (1  -0  +  P«f)  (3.18) 

In  this  case,  mean  and  fluctuating  quantities  are  fully  prescribed  by  f  and  the  governing 
equation  for  g  does  not  have  to  be  solved.  In  addition  to  p,  which  is  required  to  solve 
equation  (3.2),  time-averaged  void  fractions  are  an  important  scalar  property  to  be 
tested  using  present  measurements.  Noting  that  aQ  =  1  and  a„  =  0,  time-averaged 
mean  and  fluctuating  void  fractions  can  be  found  from  equations  (3.15)  and  (3.17)  as 
follows: 


a  =  p  f  /  p0  (3.19) 

a’2  =  p2  f  (l  -  f)  /  (p0  Poo)  (3.20) 

In  these  equations  state  o  should  be  replaced  by  state  e  if  the  effective-adapted-jet 
approximation  is  used.  Use  of  either  the  clipped  Gaussian  or  double  delta  function 
PDF's  gave  nearly  the  same  results;  thus  the  latter  is  preferred  for  simplicity  for 
noncondensing  flows. 

3.6  Computations 

Equations  (3.2)  were  integrated  using  a  modified  version  of  the  GENMIX 
algorithm  due  to  Spalding  (1977).  The  large  density  variations  of  the  present 
multiphase  flows  required  a  finer  numerical  grid  than  is  needed  for  single-phase  flows. 
Results  reported  in  the  following  used  360  cross-stream  grid  nodes,  with  streamwise 
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step  sizes  chosen  to  be  less  than  0.2  percent  of  the  current  flow  width,  or  an 
entrainment  increase  of  less  than  0.2  percent,  whichever  was  limiting.  Doubling  the 
number  of  nodes  in  the  mesh  resulted  in  less  than  a  1  percent  variation  of  flow 
properties;  therefore,  numerical  closure  was  adequate,  particularly  in  view  of  the  other 
uncertainties  of  the  analysis. 

Sensitivity  analysis  was  performed  to  find  the  influence  of  initial  conditions  on 
computed  flow  properties.  Results  for  round  jets  are  summarized  in  Table  4.  Changes 
of  100  percent  for  ke,  ee,  UCe  and  Re  have  a  negligible  effect  on  computed  results. 
Similar  changes  of  de  and  Ri  have  a  larger  effect,  however,  uncertainties  of  these 
quantities  are  much  less  than  100  percent. 


4.1 


Flow  Visualization.  Results  of  flow  visualization  for  the  round  jets  are 
discussed  by  Loth  &  Faeth  (1987)  and  Loth  (1988).  Motion  pictures  of  the  flow  (at 
roughly  1000  pictures/s)  gave  evidence  of  large-scale  unsteadiness  for  m/ms  =  0.6  and 
1.0,  similar  to  the  observations  of  others  (Kemey  et  al.,  1972;  Weimer  et  al.,  1973; 
Chan,  1974;  Cho  et  al,  1987;  Chun  &  Sonin,  1984;  Lambier  &  Chow,  1984;  Simpson 
&  Chan,  1982).  This  involved  random  fluctuations  of  gas  release  (ca.  10-20/s) 
producing  mushroom-like  gas  structures  near  the  exit  and  the  sudden  appearance  (ca.  1 
ms)  of  gas  below  the  passage  exit.  This  behavior  was  similar  to  the  pulsitile  release  of 
gas  near  the  source  of  strongly-buoyant  single-phase  plumes.  The  frequency  of 
appearance  of  large-scale  disturbances  decreased  with  increasing  m/ms  for 
underexpanded  jets.  Only  occasional  disturbances  (typically  2-3/s)  were  observed  for 
m/ms  =  8,  lasting  only  a  few  ms.  They  appear  to  be  similar  to  the  "reverse  shocks" 
mentioned  by  Surin  et  al.  (1983)  and  are  associated  with  the  external-expansion  region, 
possibly  involving  interactions  between  shock-waves  and  other  gas-dynamic 
phenomena  and  the  motion  of  gas-liquid  interfaces.  Plenum  pressure  fluctuations  due  to 
these  disturbances  decreased  monotonically  with  increasing  injector  mass  flow  rates, 
e.g.,  maximum  plenum  pressure  fluctuations  were  1.2  and  0.2  percent  at  m/ms  of  0.6 
and  8.  The  use  of  screens  at  the  passage  exit,  and  various  injector  exit  configurations, 
injector  angles  (ranging  from  horizontal  to  vertically  upward),  liquid  depths,  plenum 
sizes  and  wave-damping  methods,  had  no  observable  influence  on  the  character  of  the 
large-scale  unsteady  disturbances.  On  the  other  hand,  even  low  levels  of  liquid  coflow 
stabilized  the  disturbances  appreciably. 

Results  for  the  plane  jets  were  somewhat  similar  to  the  round  jets.  High-speed 
motion  pictures  (at  roughly  1000  pictures  per  second)  indicated  rapid  lateral  growth  of 
the  void-containing  region  of  the  flow.  Innate  unsteadiness  was  also  observed  but  to  a 
lesser  degree  than  for  the  round  jets.  Similar  to  the  round  jets,  global  unsteadiness  of 
the  plane  jets  was  reduced  as  the  underexpansion  ratio  increased.  For  adapted  flow,  the 
disturbances  were  large  enough  so  that  a  frequency  of  5.7  Hz  could  be  measured, 
however,  at  higher  underexpansion  ratios  the  disturbances  were  not  sufficiently  defined 
to  judge  their  frequency.  Since  the  exit  of  the  plane  jet  did  not  protrude  into  the  flow, 
no  reverse  flow  was  observed.  However,  preliminary  tests  with  a  protruding  plane  jet, 
having  no  sidewalls  in  the  multiphase  region,  did  produce  large  random  disturbances 
and  reverse  flow. 


Typical  photographs  of  the  appearance  of  the  plane  jet  appear  in  Fig.  7.  Figure 
7A  is  a  flash  photograph  for  air  into  water  while  Figs.  7B  and  C  are  shadowgraphs  for 
air  into  air  and  water  -  all  for  m/ms  =  3.0.  The  photograph  for  air  into  water  shows  the 
presence  of  a  series  of  oval-shaped  lobes  extending  some  distance  into  the  flow  which 
is  reminiscent  of  the  deflection  of  the  mixing  layer  for  underexpanded  gas  jets  in  gases, 
see  (Chuech  et  al.,  1988)  or  Fig.  7B.  Similarly,  comparing  Figs.  7B  and  C  indicates 
are  rather  similar  wave  structure  for  injection  of  air  into  air  and  water,  at  least  through 
the  first  lobe  (or  shock  cell).  This  provides  direct  proof  that  a  shock-wave-containing 
external-expansion  region  is  present  for  underexpanded  gas  jets  in  liquids.  The  core  of 
the  first  cell  appears  to  be  clear  of  drops  although  it  is  surrounded  by  a  multiphase 
mixing  layer.  At  the  end  of  the  first  lobe,  however,  a  drop-containing  region  is 
observed  adjacent  to  the  core  -  eminating  more  or  less  from  the  region  where  the  mixing 
layer  makes  its  closest  approach  to  the  plane  of  symmetry.  The  presence  of  drops  in  the 
core  region  tends  to  obscure  the  subsequent  lobes  but  results  in  Fig.  7A  suggest  that 
this  multiphase  gas  dynamic  region  extends  some  distance  from  the  jet  exit. 

Static  Pressures.  The  nature  of  the  external  expansion  region  for  underexpanded 
gas  jets  in  liquids  can  be  seen  more  quantitatively  from  the  static  pressure 
measurements.  Static  pressures  along  the  axis  of  round  jets,  for  m/ms  >  1.0,  are 
illustrated  in  Fig.  8.  Results  are  shown  for  injection  into  both  air  and  water,  from  the 
same  1 1.0  mm  diameter  passage.  The  results  for  injection  into  air  are  very  similar  to  the 
findings  of  Chuech  et  al.  (1988)  for  underexpanded  fully -developed  pipe  flows  in  still 
air.  Static  pressures  along  the  ais  exhibit  a  decaying  oscillatory  behavior  due  to  the 


presence  of  shock  cells  that  eventually  decay  away  as  the  mixing  layers  near  the  edge  of 
the  flow  reach  the  axis.  The  static  pressure  variations  along  the  axis  for  injection  of  air 
into  water  are  similar  to  those  observed  for  injection  of  air  into  air,  at  least  for  the  first 
few  shock  cells.  Notably,  Surin  et  al.  (1983)  observed  similar  agreement  of  oscillatory 
dynamic  pressure  variations  along  the  axis  of  underexpanded  gas  jets  in  gases  and 
liquids.  This  confirms  that  a  shock-wave-containing  external-expansion  region  is 
present  for  underexpanded  gas  jets  in  liquids.  The  main  difference  between  the  static 
pressure  records  for  underexpanded  air  jets  in  air  and  in  water  is  that  the  external- 
expansion  region  d  •'-ays  more  rapidly  for  injection  into  water.  This  behavior  is 
expected  since  turbulent  mixing  is  more  rapid  when  a  low-density  material  is  injected 
into  a  high-density  environment  (Chen  &  Faeth,  1982,  1983;  Faeth,  1987). 

Static  pressures  along  the  axis  of  the  plane  jets  are  illustrated  in  Fig.  9  for  m/ms 
>  1.0  for  injection  into  both  air  and  water.  The  results  are  generally  similar  to  the 
findings  for  the  round  jets:  the  static  pressure  variation  in  the  external-expansion  region 
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Plane  injector  at  m/mg  =  3.0:  A)  flash  photograph  of  air  into 
water,  B)  shadowgraph  of  air  into  air,  C)  shadowgraph  of  air 
into  water. 
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Figure  8.  Static  pressure  distribution  along  axis  for  round  jets. 
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is  similar  for  injection  into  air  and  water  although  the  latter  decays  more  rapidly. 
Computations  using  a  parabolized  Navier-Stokes  method  described  by  Chuech  (1987) 
and  Cheuch  et  al.  (1988)  are  also  shown  on  the  plot  for  comparison  with  measured 
properties  of  air  jets  in  air.  The  comparison  between  predictions  and  measurements  is 
not  very  good,  largely  due  to  the  presence  of  acoustic  feedback  which  is  discussed 
next. 

Acoustic  Feedback.  Another  important  feature  of  the  passage  exit  condition  is 
the  presence  of  acoustic  feedback  within  the  shock-cell  pattern.  Acoustic  feedback 
arises  as  a  result  of  pressure  disturbances  traveling  upstream  near  the  edge  of  the  jet, 
reflecting  from  the  surfaces  near  the  exit  plane  of  the  nozzle  and  resonating  at  a 
particular  frequency.  The  tuning  condition  is  set  by  the  sound  wave  speed  of  the 
ambient  fluid  and  an  acoustic  reflection  distance  from  the  jet  exit  to  some  downstream 
position  of  the  external  expansion  region,  associated  with  one  of  the  shock  cells. 
Sherman  et  al.  (1976)  have  reported  substantially  increased  mixing  due  to  acoustic 
feedback  for  underexpanded  air  jets  in  air,  with  dominant  feedback  frequencies 
approximately  proportional  to  the  sonic  speed  divided  by  the  first  or  second  shock-cell 
length. 

Measurements  with  a  hydrophone  at  a  radial/lateral  position  of  six  exit  widths 
from  the  centerline  and  six  exit  widths  downstream  were  completed  for  the  4.9  mm 
round  injector  and  the  4.8  mm  plane  injector,  for  air  injected  both  into  air  and  into 
water.  Power  spectral  densities  of  the  acoustic  signal  are  plotted  as  a  function  of 
frequency  in  Figs.  10  and  1 1  for  the  round  and  plane  jets,  respectively.  Measurements 
were  made  for  frequencies  of  up  to  62  kHz  (the  approximate  Strouhal  number),  with 
significant  signal  content  found  in  the  1-25  kHz  regime.  Little  evidence  of  acoustic 
feedback  for  the  round  jet  into  still  air  was  observed,  except  for  a  resonance  at  23  kHz 
for  a  mass  flow  ratio  of  two.  This  is  not  surprising  since  the  injector  was  protruding 
and  provided  little  surface  area  to  reflect  acoustic  waves.  For  injection  into  still  water, 
signal  levels  are  typically  lower  at  high  frequencies,  ca.  5-25  kHz,  and  no  evidence  of 
acoustic  feedback  is  observed.  This  behavior  may  be  due  to  the  higher  acoustic  velocity 
of  water  than  air,  which  would  vastly  modify  tuning  conditions  for  acoustic  feedback. 

There  is  substantial  evidence  of  acoustic  feedback  in  Fig.  11  for  all  the 
underexpansion  ratios  for  the  plane  air  jet  into  air.  Feedback  peaks  are  especially  large 
at  mass  flow  ratios  of  two  and  three,  with  frequencies  corresponding  to  acoustic 
reflection  distances  of  several  exit  widths.  The  flat  surface  along  the  exit  plane  is 
probably  responsible  for  this  dramatic  increase  in  acoustic  feedback,  since  it  provides  a 
reflecting  surface  for  acoustic  disturbances.  These  strong  effects  of  acoustic  feedback 
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probably  explain  the  more  rapid  decay  of  static  pressures  than  predicted  in  Fig.  9  for 
underexpanded  air  jets  in  air.  Similar  to  the  round  jet,  however,  injection  into  water 
yields  lower  acoustic  levels  at  the  high  frequencies,  ca.  5-25  kHz,  and  no  evidence  of 
acoustic  feedback.  In  addition  to  effects  of  the  different  acoustic  velocities  of  water  and 
air,  the  presence  of  bubbles  in  the  multiphase  flow  may  also  damp  pressure  waves,  as 
reported  by  Borisov  (1983).  These  acoustic  measurements  indicate  that  the  plane  jet 
may  exhibit  enhanced  mixing  in  comparison  to  the  round  injector  due  to  the  presence  of 
a  reflective  surface  near  the  jet  exit,  similar  to  the  findings  of  Sherman  et  al.  (1976). 
Present  measurements  also  indicate  that  acoustic  feedback  is  eliminated  when  water  is 
the  ambient  medium,  even  when  reflective  surfaces  are  present 

4.2  Flow  Structure 
4.2.1  Round  Jets 

Initial  Conditions.  Measurements  of  mean  and  fluctuating  velocities  at  the  exit 
of  the  round  jets,  for  all  underexpansion  ratios  considered  here,  are  reported  by  Loth  & 
Faeth  (1987)  As  noted  earlier,  these  results  confirmed  that  the  flows  were  essentially 
fully-developed  pipe  flows  at  the  jet  exit. 

Void  Fractions.  Measurements  and  predictions  using  standard  initial  conditions 
for  round  jets  have  already  been  discussed  by  Loth  &  Faeth  (1987).  These  results  will 
be  revisited  in  the  following,  considering  the  effect  of  the  using  of  mixed  initial 
conditions  and  well  as  additional  measurements  of  Bell  et  al.  (1972),  Carreau  et  al. 
(!985)  and  Tross  (1974). 

Time-averaged  void  fractions  along  the  axis  of  the  round  jets,  otc  are  illustrated 
in  Fig.  12.  Results  shown  include  present  measurements  for  d  =  4.9  mm,  the  probe 
measurements  of  Tross  (1974)  for  m/ms  =1.0  and  2.0,  and  predictions  (denoted  EEC 
model)  using  both  standard  and  mixed  initial  conditions.  The  measurements  of  Tross 
(1974)  are  consistently  lower  than  the  present  measurements  in  the  region  where  they 
can  be  compared.  This  behavior  is  expected  since  all  biases  of  probes  tend  to  reduce 
observed  void  fractions  (Tross,  1974).  It  is  felt  that  the  present  nonintrusive  gamma- 
ray  absorption  measurements  are  more  reliable.  The  measurements  show  that  mixing 
rates  along  the  axis  are  progressively  reduced  as  the  underexpansion  ratio  is  increased. 
This  slower  rate  of  mixing  is  the  price  that  must  be  paid  to  obtain  more  stable  injector 
operation  by  increasing  m/«ns. 


Predictions  agree  reasonably  well  with  present  measurements  in  Fig.  12,  for 
m/ms  >  2  where  effects  of  unsteadiness  are  reduced.  Use  of  mixed  initial  conditions 


provide  somewhat  better  agreement  with  measurements  than  the  standard  initial 
conditions  at  lower  values  of  m/ms,  however,  both  methods  still  underestimate  mixing 
rates  at  m/ms  =  0.6.  This  deficiency  cannot  be  due  to  the  LHF  approximation,  which 
invariably  causes  mixing  rates  of  multiphase  flows  to  be  overestimated  (Faeth,  1987) 
while  the  effective-adapted-jet  approximation  is  not  used  at  m/ms  =  0.6.  Thus, 
increased  large-scale  unsteadiness  observed  at  low  m/ms  is  probably  responsible  foi 
enhanced  mixing  exhibited  by  the  measurements. 

Measured  and  predicted  radial  distributions  of  time-averaged  void  fractions,  a, 
are  illustrated  in  Figs.  13-17  for  mJms  =  0.6,  1.0,  2.0,  4.0  and  8.0.  The  results  are 
plotted  as  a  function  of  r/x,  which  is  the  radial  similarity  variable  for  fully-developed 
single-phase  turbulent  jets,  so  that  the  actual  width  of  the  flow  can  be  seen.  In  addition 
to  the  present  data  for  4.9  and  1 1 .0  mm  diameter  passages,  probe  measurements  from 
Tross  (1974)  are  also  shown  for  m/ms  =1.0  and  2.0.  Similar  to  the  measurements 
along  the  axis,  the  measurements  of  Tross  (1974)  are  generally  lower  than  present 
measurements,  probably  due  to  probe  biases. 

The  most  striking  feature  of  the  measurements  illustrated  in  Figs.  13-17  is  the 
unusual  width  of  the  void  fraction  profiles.  This  is  expected  near  the  passage  exit  when 
results  are  plotted  as  a  function  of  r/x.  Furthermore,  at  x/d  =  4  and  m/ms  =  2.0  and  4.0, 
there  is  a  bulge  in  the  void  fraction  distribution  near  the  edge  of  the  flow  which  is 
probably  associated  with  the  presence  of  the  external-expansion  region  at  this  position 
(see  Fig.  8).  Nevertheless,  flow  widths  remain  unusually  large  even  at  x/d  =  80  where 
the  flow  might  be  expected  to  approach  the  properties  of  fully-developed  jets,  e.g.,  the 
edge  of  the  gas-containing  region  is  at  r/x  ~  0.3  at  x/d  =  80,  which  is  nearly  twice  the 
width  observed  for  properties  in  fully-developed  single-phase  jets  (Wygmanski  & 
Fiedler,  1969).  The  main  reason  for  this  behavior  is  that  void  fraction  is  an  unusually 
sensitive  indicator  of  the  mixing  level  of  the  flow  at  low  void  fractions  (see  Fig.  6). 
Recall  that  when  the  mixture  fraction  f  <  <  pa/pw»  the  density  ratio  of  air  and  water,  the 
relationship  between  void  fraction  and  mixture  fraction  becomes  (equation  (3.6)): 

a  =  pwf/pa  (4.1) 

For  present  test  conditions  Pw  /  Pa  ~  850;  therefore,  a  is  still  significant  even  when  f  is 
much  smaller  than  values  of  mixing  levels  that  would  normally  be  associated  with  the 
edge  of  a  jet. 

In  view  of  the  sensitivity  of  void  fractions  to  mixing  levels  (mixture  fraction), 
the  comparison  between  predictions  and  measurements  in  Fig.  13-17  is  encouraging,  in 
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Figure  13.  Time-averaged  void  fractions  of  round  jets  for  m/me 
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Figure  14.  Time-averaged  void  fractions  of  round  jets  for  m/ms 
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Figure  16.  Time-averaged  void  fractions  of  round  jets  for  m/ms  -  4.0. 
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the  sense  that  predicted  characteristic  flow  widths  tend  to  scale  with  the  measurements 
as  x/d  and  m/ms  vary.  However,  the  evolution  of  the  profiles  from  the  distribution 
prescribed  at  the  jet  exit  is  relatively  slow  and  profiles  using  standard  initial  conditions 
are  far  too  blunt  in  the  radial  direction,  particularly  near  the  passage  exit.  Use  of  the 
mixed  initial  conditions  yields  much  better  predictions  of  the  shape  of  the  profiles, 
particularly  near  the  jet  exit.  This  suggests  that  mixing  is  enhanced  near  the  jet  exit  for 
all  conditions,  probably  as  a  result  of  intrinsic  unsteadiness.  In  particular,  failure  of  the 
LHF  approximation  generally  causes  overestimation  of  flow  widths  (Faeth.  1987), 
thus,  separated- flow  phenomena  do  not  provide  an  obvious  explanation  of  the 
discrepancies  between  predictions  and  measurements  seen  in  Figs.  13-17.  Similarly, 
the  shapes  of  void  fraction  profiles  do  not  differ  for  adapted  and  underexpanded  jets; 
thus,  use  of  the  effective-adapted-jet  approximation  is  not  an  obvious  source  of  the 
discrepancies  either. 

Figure  18  is  an  illustration  of  measured  and  predicted  characteristic  void- 
fraction-diameters,  2ro.5  (defined  as  twice  the  radius  where  oc/otc  =  0.5)  plotted  as  a 
function  of  x/d  for  various  m/ms.  Also  shown  on  the  plot  is  the  radius  of  the  gas 
containing  region,  which  is  roughly  comparable  to  2ro.5,  measured  by  Bell  et  al.  (1972) 
using  a  passage  plate  having  various  hole  diameters  for  an  underexpansion  ratio  of  3.7. 
The  results  of  Bell  et  al.  (1972)  are  in  reasonably  good  agreement  with  present 
measurements  in  view  of  the  somewhat  different  definitions  of  characteristic  diameters. 
The  unusual  width  of  these  flows,  when  viewed  in  terms  of  a,  is  evident  from  the 
figure.  As  before,  this  reflects  the  sensitivity  of  void  fractions  to  mixing  levels  (note 
that  a  void  fraction  of  0.5  corresponds  to  mixture  fractions  of  roughly  0.1  percent).  The 
characteristic  void  fraction  diameters  increase  everywhere  as  the  underexpansion  ratio 
increases:  even  at  x/d  =  80,  2ro.5  /  x  ~  0.2  for  m/ms  =  0.6  and  ~  0.3  for  m/ms  =  8.0. 
Predictions  using  the  two  initial  conditions  only  differ  near  the  jet  exit  and  are  in 
reasonably  good  agreement  with  the  measurements. 

Dynamic  Pressures.  Other  flow  properties  exhibit  more  conventional 
characteristic  flow  widths  than  void  fractions  since  they  are  less  sensitive  to  mixing 
levels.  This  can  be  illustrated  using  measurements  of  dynamic  pressures  from  Tross 
(1974)  for  conditions  very  similar  to  the  present  investigation.  Although  these 
measurements  were  obtained  using  probes,  reasonable  comparisons  between  the  probe 
measurements  and  present  nonintrusive  measurements  of  void  fractions  provide  some 
confidence  in  the  results. 

Measured  and  predicted  dynamic  pressures  along  the  axis  are  illustrated  in  Fig. 
19.  The  results  involve  two  initial  jet  diameters,  3.2  and  6.4  mm,  and  m/ms  =  1.0,  1.5 
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and  2.0.  Although  void  fractions  do  not  increase  appreciably  from  unity  until  x/d  ~  30 
for  these  conditions  (see  Fig.  12),  dynamic  pressures  begin  to  decrease  before  x/d  =  10 
and  then  exhibit  a  logarithmic  decay  typical  of  fully-developed  single-phase  jets.  This 
difference  is  a  result  of  the  large  density  ratio  of  the  flow,  i.e.,  small  volumes  of  liquid 
can  absorb  large  amounts  of  momentum  from  the  gas.  Predictions  based  on  the 
standard  initial  conditions  tend  to  overestimate  the  rate  of  decay  of  dynamic  pressures; 
in  contrast,  use  of  mixed  initial  conditions  yields  reasonably  good  agreement  with  the 
measurements  -  supporting  the  presence  of  unusually  high  mixing  levels  near  the  jet 
exit.  An  obvious  defect  of  the  predictions,  however,  is  that  use  of  the  effective- adapted- 
jet  approximation  cannot  provide  information  concerning  dynamic  pressure  variations 
in  the  external-expansion  region,  analogous  to  the  static  pressures  illustrated  in  Fig.  8. 

Radial  distributions  of  dynamic  pressures  reported  by  Tross  (1974)  are 
illustrated  in  Figs.  20-22  for  m/ms  =  1.0,  1.5  and  2.0.  Predictions  using  standard  and 
mixed  initial  conditions  are  also  shown  on  the  plots:  there  is  little  difference  between  the 
two  predictions  for  these  variables. 

Flow  widths  exhibited  by  dynamic  pressures  in  Figs.  20-22  are  far  more  typical 
of  the  properties  of  single-phase  jets  than  flow  widths  based  on  void  fractions,  e.g.,  the 
edge  of  the  flow  is  at  r/x  ~  0.15  for  x/d  =  72.  Present  predictions  and  the  measurements 
are  in  good  agreement  for  the  underexpanded  jets:  m/ms  =  1.5  and  2.0.  Predictions 
indicate  a  narrower  flow  field  than  measured  for  the  adapted  jet  (Fig.  20).  This  is 
consistent  with  the  presence  of  significant  effects  of  unsteadiness  at  this  condition, 
which  would  be  expected  to  enhance  mixing  rates. 

Entrainment.  Predicted  and  measured  entrainment  coefficients  for  the  round 
jets,  defined  according  to  equation  (2.2)  are  plotted  in  Fig.  23.  Experimental  results 
include  present  measurements:  measurements  of  Bell  et  al.  (1972)  and  Carreau  et  al. 
(1985)  for  horizontal  underexpanded  nitrogen  jets  in  water,  using  the  water  overflow 
technique;  and  the  empirical  correlation  of  Ricou  &  Spalding  (1961)  for  the  fully- 
developed  region  of  single-phase,  variable-density  gas  jets.  The  present  measurements 
and  those  of  Carreau  et  al.  (1985)  are  in  excellent  agreement  The  earlier  measurements 
of  Bell  et  al.  (1972),  however,  are  more  scattered  and  are  generally  higher  than  the  rest: 
the  reason  for  this  difference  in  behavior  is  not  known.  Measured  entrainment 
coefficients  are  relatively  constant  at  small  x/d  but  begin  to  increase  at  larger  x/d,  with 
the  increase  beginning  sooner  as  m/ms  is  reduced.  This  is  an  effect  of  buoyancy 
increasing  the  momentum  of  the  flow,  which  becomes  important  nearer  to  the  jet  exit  as 
m/ms  decreases,  since  this  reduces  the  initial  streamwise  momentum  of  the  flow.  When 
effects  of  buoyancy  are  small,  the  entrainment  coefficient  remains  nearly  constant  in  the 
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Figure  20.  Mean  dynamic  pressures  of  round  jets  for  m/ms  =  1.0. 
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fully-developed  portion  of  the  flow,  similar  to  the  correlation  of  Ricou  &  Spalding 
(1961).  Measured  entrainment  coefficients  also  increase  everywhere  as  m/ms  decreases. 
This  probably  due  to  enhanced  mixing  caused  by  increasing  unsteadiness  at  low 
underexpansion  ratios. 

The  comparison  between  predictions  using  standard  initial  conditions  and  the 
measurements,  seen  in  Fig.  23,  is  generally  poor,  although  predictions  tend  to  improve 
as  m/ms  and  x/d  increase.  Predictions  using  the  mixed  initial  conditions,  however,  are 
in  reasonably  good  agreement  with  the  measurements  except  for  m/ms  =  0.6  where 
predictions  substantially  underestimate  the  measurements.  The  behavior  suggests  that 
the  discrepancies  are  largely  due  to  enhanced  mixing  from  effects  of  unsteadiness  very 
near  the  jet  exit.  This  is  a  major  effect  when  the  flow  at  the  jet  exit  is  subsonic  but  still 
persists  to  a  significant  degree  for  the  underexpanded  jets. 


4.2.2  Plane  Jets 

Initial  Conditions.  Time  averaged  mean  and  fluctuating  axial  velocities  were 
also  measured  for  the  plane  jet  in  still  air.  The  results  are  illustrated  in  Figs.  24  and  25. 
The  plane  converging  nozzle  was  designed  to  yield  a  uniform  exit  velocity  profile.  The 
measured  mean  velocity  profile  is  reasonably  uniform  except  for  the  last  ten  percent 
near  the  jet  edge.  The  variation  of  mean  velocities  near  the  jet  edge  is  probably  due  to 
gradient  broadening  of  the  LDA  measurements,  since  the  size  of  the  measuring  volume 
is  approximately  equal  to  the  distance  between  data  points  in  this  region.  Note  that  since 
measurements  were  taken  at  an  x/b  position  of  zero,  the  expansion  fan  acceleration 
observed  for  the  round  jet  (Loth  &  Faeth,  1987)  is  not  observed.  Results  illustrated  in 
Fig.  25  indicate  a  uniform  turbulence  level  of  about  three  percent  but  the  spatial 
resolution  is  not  adequate  to  detect  the  edge  boundary  layer.  From  these  measurements, 
the  nozzle  appears  to  provide  a  uniform  velocity  profile  with  a  relatively  low  turbulence 
intensity  at  the  jet  exit 

Void  Fractions.  Measurements  of  void  fraction  distributions  were  limited  to  the 
region  1  £  x/b  <  8.  Measurements  were  not  undertaken  at  larger  x/b  since  the  rapid 
growth  of  the  flow  caused  the  flow  aspect  ratio  to  decline  to  two  at  x/b  =  8  so  that  two- 
dimensionality  is  questionable  farther  downstream.  Measurements  and  all  predictions 
indicated  a  time-averaged  void  fraction  along  the  axis  of  unity  in  the  region  where 
measurements  were  made;  therefore,  only  void  fraction  distributions  at  various  x/b  will 
be  considered  in  the  following. 

Measured  and  predicted  distributions  of  time-averaged  void  fractions  for  m/ms 
=  1,  2,  3  and  4  are  illustrated  in  Figs.  26-29.  Results  are  plotted  as  a  function  of  y/x, 
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Figure  27.  Time-averaged  void  fractions  of  plane  jets  for  m/ms 
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Figure  28.  Time-averaged  void  fractions  of  plane  jets  for  m/ms  =  3. 
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Figure  29,  Time-averaged  void  fractions  of  plane  jets  for  m/ms  -  4. 


which  is  the  similarity  variable  for  the  fully-developed  region  of  plane  turbulent  jets  and 
plumes.  The  predictions  are  based  on  the  equivalent-adapted-jet  approach  (EEC  model) 
using  both  standard  and  mixed  initial  conditions. 

Similar  to  the  round  jet  results,  the  most  striking  feature  of  the  results  illustrated 
in  Figs.  26-29  is  the  large  lateral  widths  of  the  void  fraction  profiles,  ca.  y/x  of  0.5  as 
opposed  to  0.15  for  scalar  property  widths  of  single-phase  flows.  This  behavior  is 
partly  caused  by  the  fact  that  the  near  jet  exit  region  is  being  considered  where  flow 
widths  are  invariably  quite  large  when  plotted  in  terms  of  y/x.  However,  the  large 
sensitivity  of  void  fractions  to  mixture  fractions  is  also  a  factor,  as  pointed  out  in 
connection  with  the  round  jets.  In  general  flow  widths  become  narrower,  in  terms  of 
y/x,  with  increasing  x/b.  This  follows  since  the  flow  is  still  developing  in  the  near 
injector  region,  e.g.,  the  round  jet  results  in  Figs.  13-18  suggest  that  similarity  would 
be  approached  at  larger  x/b. 

The  comparison  between  predictions  and  measurements  for  the  sonic  adapted 
jet,  Fig.  26,  is  quite  encouraging.  Use  of  both  the  standard  and  mixed  initial  conditions 
yield  a  reasonably  good  estimate  of  the  mean  void  fraction  width  of  the  flow  and  its 
development  with  increasing  distance  from  the  jet  exit.  Use  of  the  mixed  initial 
conditions,  however,  provides  a  superior  estimate  of  the  shape  of  the  void  fraction 
profile.  This  suggests  enhanced  mixing  very  near  the  injector  exit,  related  to  the  near 
injector  unsteadiness  that  was  observed  for  both  the  round  and  plane  jets.  Once  this  is 
accounted  for,  however,  subsequent  mixing  rates  are  not  unusual  in  comparison  to 
estimates  based  on  the  locally  homogeneous  flow  approximation.  This  suggests  that 
finite-rate  interphase  transport  rates  are  not  a  major  factor  with  respect  to  mixing  rates  in 
the  near-injector  region  for  adapted  conditions. 

Void  fraction  distributions  for  underexpanded  conditions,  shown  in  Figs.  27- 
29,  involve  additional  complications  due  to  the  presence  of  the  external  expansion 
region.  This  introduces  new  sources  of  errors  between  predictions  and  measurements 
due  to  varying  pressures  in  the  flow  field,  seen  in  Fig.  9,  which  cause  crosstream 
expansion  and  contraction  of  the  mean  position  of  the  mixing  layer.  Naturally  such 
features  cannot  be  approximated  by  the  adapted  jet  approximation  even  though  this 
approach  provides  a  crude  representation  of  the  mean  position  of  the  mixing  layer. 

Effects  of  the  wave  field  of  the  external  expansion  region  on  crosstream  mean 
void  fraction  distributions  are  evidenced  by  the  apparently  whimsical  variations  in  the 
agreement  between  predictions  and  measurements  (particularly  when  using  the  mixed 
initial  conditions)  as  a  function  of  distance  from  the  jet  exit  and  underexpansion  ratio. 
For  example,  the  agreement  between  predictions  and  measurements  is  generally  quite 
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good  for  all  underexpansion  ratios  at  x/b  =  2;  predictions  generally  overestimate 
measured  flow  widths  at  x/b  =  1;  and  predictions  generally  underestimate  measured 
flow  widths  at  x/b  =  4  and  8.  Furthermore,  discrepancies  at  x/b  =  1,4  and  8  tend  to 
increase  as  the  underexpansion  ratio  increases.  The  reason  for  this  behavior  can  be  seen 
from  the  static  pressure  traces  illustrated  in  Fig.  9.  It  happens  that  x/b  =  1,4  and  8  all 
fall  near  conditions  where  static  pressure  near  the  plane  of  symmetry  is  near  the  ambient 
pressure  and  is  decreasing  with  increasing  distance  from  the  jet  exit.  This  corresponds 
to  positions  where  the  outward  deflection  of  the  mixing  layer  surrounding  the  external 
expansion  region  is  a  maximum  causing  the  flow  to  be  wider  than  predictions  ignoring 
the  pressure  field.  Naturally  this  outward  deflection  increases  as  the  underexpansion 
ratio  increases,  so  that  discrepancies  between  predictions  and  measurements  increase 
with  increasing  underexpansion  ratio.  On  the  other  hand,  the  position  x/b  =  2 
corresponds  to  a  minimum  in  the  static  pressure  field  along  the  axis,  where  the  mixing 
layer  is  roughly  halfway  between  its  extreme  positions  -  which  fortuitously 
corresponds  to  the  location  predicted  when  the  pressure  field  is  ignored. 

Lateral  deflections  of  the  mixing  layers  due  to  pressure  variations  can  only  be 
predicted  by  extending  the  analysis  to  consider  the  wave  field,  analogous  to  parabolized 
Navier-Stokes  methods  for  underexpanded  air  jets  in  air  (Chuech  et  al.,  1988). 
However,  use  of  the  locally  homogeneous  flow  approximation  introduces  very  large 
density  variations  in  the  flow  (ca.  1000:1)  which  makes  it  much  more  difficult  to 
achieve  stable  and  accurate  solutions  of  the  governing  equations.  Since  parabolized 
Navier  Stokes  algorithms  also  require  careful  handling  to  assure  stable  and  accurate 
results  for  even  single-phase  flows,  extension  of  these  methods  to  treat  gas  injection 
into  liquids  under  the  locally  homogeneous  flow  approximation  will  not  be  a  trivial 
matter  -  difficulties  concerning  the  proper  treatment  of  compressibility  effects  and  shock 
wave  interactions  in  a  multiphase  flow  aside. 

In  addition  to  discrepancies  just  noted  due  to  the  external  expansion  region, 
differences  of  predictions  between  the  standard  and  mixed  initial  conditions  are  roughly 
the  same  as  for  the  round  jets,  with  use  of  the  mixed  initial  conditions  providing 
superior  predictions  of  the  general  shape  of  the  void  fraction  mixing  layer.  This  adds  to 
evidence  for  unusual  degrees  of  mixing  near  the  jet  exit,  although  the  mechanism  for 
this  behavior  has  not  yet  been  resolved. 

Figure  30  is  an  illustration  of  predicted  and  measured  characteristic  flow  widths 
(defined  as  twice  the  crosstream  distance  where  a/cxc  =  0.5)  plotted  as  a  function  of  x/b 
for  various  underexpansion  ratios.  These  results  illustrate  the  large  width  of  the  flow 
near  the  jet  exit,  which  is  expected  due  to  flow  development,  with  a  trend  toward  an 


O  DATA _ 

EEC  MODEL _ 

-  STANDARD  I.  C. 

- MIXED  I.C. 


asymptotic  condition  where  the  flow  width  is  relatively  independent  of  the  distance 
from  the  jet  exit  at  larger  x/b.  As  noted  earlier,  however,  this  far  field  flow  width  is 
unusually  large  due  to  the  influence  of  the  large  density  ratio  of  the  flow  on  the 
sensitivity  of  the  void  fractions  to  mixing  levels.  Similar  to  results  for  the  round  jets, 
the  present  approximate  analysis  yields  reasonably  good  estimates  of  the  characteristic 
flow  widths  and  trends  with  respect  to  distance  from  the  jet  exit  and  the  underexpansion 
ratio. 

Entrainment.  Predicted  and  measured  dimensionless  entrainment  coefficients, 
based  on  equation  (2.2),  are  plotted  as  a  function  of  x/b  for  various  underexpansion 
ratios  in  Fig.  31.  An  emperical  relationship  for  gas  injection  into  a  still  gas,  due  to 
Schneider  (1985),  is  also  plotted  on  the  figure.  Measurements  exhibit  a  slight  reduction 
in  entrainment  rates  for  increasing  underexpansion  ratios.  This  may  be  due  to  increased 
stability  of  the  near  injector  region  to  flow  unsteadiness,  which  was  also  noted  in 
connection  with  the  measurements  for  underexpanded  round  jets,  see  Fig.  23. 

Another  factor  is  that  Mach  numbers  in  the  external  expansion  region  increase 
as  the  underexpansion  ratio  increases:  this  would  tend  to  reduce  mixing  rates  analagous 
to  well  known  effects  of  compressibility  seen  in  supersonic  shear  layers  (Papamoschou 
&  Rosko,  1986;  Bogdanov,  1983).  The  measurements  also  exhibit  a  progressive 
reduction  in  the  dimensionless  entrainment  rate  with  increasing  x/b,  at  all 
underexpansion  ratios.  This  effect  is  largely  due  to  transition  of  flow  properties  from 
nearly  slug  flow  at  the  jet  exit  to  fully-developed  flow  at  downstream  locations  and  can 
be  found  from  the  scaling  laws  for  a  fully-developed  turbulent  jet.  Indeed,  the  empirical 
relationship  due  to  Schneider  (1985)  for  gas  jets  in  still  gases  exhibits  a  similar  trend, 
although  quantitative  agreement  between  this  relationship  and  the  measurements  is  not 
very  good. 

The  comparison  between  predicted  and  measured  dimensionless  entrainment 
rates  in  Fig.  31,  using  the  standard  initial  conditions,  is  generally  not  very  good.  This 
behavior  can  probably  be  attributed  to  unsteady  features  of  the  flow  near  the  jet  exit. 
This  interpretation  is  supported  by  the  fact  that  predictions  using  the  mixed  initial 
conditions  provide  very  encouraging  agreement  with  measurements,  including  trends 
with  respect  to  both  underexpansion  ratio  and  distance  from  the  jet  exit.  This  result 
emphasizes  the  sensitivity  of  the  flow  to  initial  conditions  and  disturbances  near  the  jet 
exit. 


,ion 


Major  conclusions  of  the  present  study  are  as  follows: 

Shock-wave-containing  external  expansion  regions  are  present  for 
underexpanded  air  jets  in  water  but  are  smaller  in  extent  than  air  jets  in  air  for 
comparable  conditions  due  to  faster  mixing  rates.  This  behavior  was  established 
by  both  static  pressure  measurements  and  shadowgraphs  which  showed  wave 
structure,  complimenting  earlier  findings  by  Moiseev  (1962)  and  Surin  et  al. 
(1983)  based  on  dynamic  pressure  measurements. 

Flow  visualization  showed  that  the  mixing  layer  surrounding  the  external 
expansion  region  for  plane  air  jets  in  water  involved  a  bubbly  region  toward  the 
liquid  side  and  a  drop-containing  region  toward  the  gas  side,  although 
significant  stripping  of  drops  only  began  near  the  downstream  end  of  the  first 
shock  cell. 

Effects  of  acoustic  feedback  observed  for  underexpanded  air  jets  in  air 
disappeared  for  underexpanded  air  jets  in  water  at  comparable  conditions.  This 
behavior  is  probably  due  to  changed  shock  wave  positions  in  the  external- 
expansion  region  and  different  acoustic  velocities  and  acoustic  dampening  in  the 
multiphase  mixing  layer. 

Air  jets  in  water  exhibit  unusual  flow  widths  based  on  void  fractions  due  to  the 
strong  sensitivity  of  void  fractions  to  mixing  levels  (mixture  fraction),  however, 
flow  diameters  based  on  dynamic  pressures  in  round  jets  are  more 
conventional. 

Use  of  the  locally-homogeneous-flow  and  effective-adapted-jet  approximations 
yielded  encouraging  predictions  of  the  properties  of  round  air  jets  in  water, 
particularly  at  large  m/ms  where  effects  of  unsteadiness  were  reduced.  The 
prescription  of  initial  conditions  for  these  calculations,  however,  had  a 
significant  influence  on  predictions  of  the  shape  of  radial  void  fraction 
distributions  and  entrainment  coefficients.  Use  of  mixed  initial  conditions,  to 
account  for  enhanced  mixing  due  to  unsteadiness  near  the  jet  exit,  yielded  best 
agreement  with  measurements. 

Use  of  the  locally-homogeneous-flow  and  effective-adapted-jet  approximations 
were  somewhat  less  successful  for  plane  jets,  where  current  measurements 
placed  greater  emphasis  on  the  near-jet-exit  region  (x/b  <  8).  The  main 
difference  is  that  predicted  lateral  profiles  of  time-averaged  void  fractions  were 
less  satisfactory  due  to  flow  deflection  by  shock  cells  in  the  external  expansion 
region  which  cannot  be  treated  using  the  effective  adapted  jet  approximation. 


While  the  structure  and  mixing  properties  of  the  vertical  underexpanded  jets 
considered  during  the  present  study  and  by  Tross  (1974),  and  the  horizontal 
underexpanded  jets  considered  by  Carreau  et  al.  (1985),  are  predicted  reasonably  well 
using  the  mixed  initial  conditions,  use  of  this  prescription  to  allow  for  near-injector 
unsteadiness  should  be  approached  with  caution.  In  particular,  Surin  et  al.  (1983)  point 
out  that  condensing  jets  exhibit  less  unsteadiness  than  noncondensing  jets  while  Loth 
(1988)  found  that  even  low  levels  of  liquid  coflow  can  significantly  stabilize  the  flow  as 
well.  Clearly,  a  better  understanding  of  the  mechanism  of  unsteadiness  near  the  jet  exit, 
in  the  presence  of  a  shock-wave-containing  external-expansion  region,  is  required  to 
properly  understand  and  reliably  predict  the  properties  of  these  flows. 
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